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Abstract 

Scattering properties of a material are changed when the material is injected 
with small acoustically soft particles. It is shown that its new scattering behavior 
can be understood as a solution of a potential scattering problem with the potential 
q explicitly related to the density of the small particles. In this paper we examine 
the inverse problem of designing a material with the desired focusing properties. 
An algorithm for such a problem is examined from the theoretical as well as from 
the numerical perspective. 



1 Introduction 

Let D C M 3 be a bounded connected domain with Lipschitz boundary S. Denote by 
no(x) the refraction coefficient in D, uq{x) = 1 in D' := M 3 \ D. Then the scattering of 
a plane acoustic wave uq = uq(x) = e tka ' x , incident upon D, is described by the system: 

[V 2 + k 2 n (x)]u{x) = in R 3 , (1) 
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u(x) = uq(x) + v(x), 



(2) 



e ikr / y \ x 

v(x) = A(a', a) ho-, r := \x\ — > oo, — := o/, (3) 

r \r J r 

where t> (x) is the scattered field, a <E S 2 is the direction of the incident plane wave, and 
a' is the direction of the scattered wave. The coefficient A(a', a) is called the scattering 
amplitude, k > is the wave number, which is assumed to be fixed throughout the paper. 
For this reason the dependence of A on k is not shown. 
Let D m , 1 < m < M, be a small particle, i.e., 

fc fl C 1, where a = - max diamD m , /c = fcmax |n (a;)|. (4) 

2 l<m<M x£D 

The geometrical shape of D m is arbitrary, but we assume that each D m has a Lipschitz 
boundary. Moreover, the Lipschitz constant is the same for every domain D m . This is 
a technical assumption which can be relaxed. It allows one to use the properties of the 
electrostatic potentials. Let 

d := min dist(D m , Dj). (5) 

Assume that 

a < d. (6) 

We do not assume that d 3> A , that is, that the distance between the particles is much 
larger than the wavelength. Under our assumptions it is possible that there are many 
small particles within the distances of the order of magnitude of the wavelength. 
The particles are assumed to be acoustically soft, i.e., 

u\ Sm = 1 < m < M. (7) 

As a result of the distribution of many small particles in D, one obtains a new material. 
We would like this "smart" material to have some desired properties. Specifically, we 
want this material to scatter the incident plane wave according to an a priori given desired 
radiation pattern, for example to focus the incident wave within a given solid angle. Is 
this possible? If yes, then how does one distribute the small particles in order to create 
such a material? In mathematical terms the problem is 

Given an arbitrary function f(f3) G L 2 (S 2 ), can one distribute small particles in D 
so that the resulting medium generates the radiation pattern A(f3) := A(/3,a), at a fixed 
k > and a fixed a <E S 2 , such that 

\\f(f3)-A(f3)\\L*(s*)<e, (8) 

where e > is an arbitrary small fixed number? 

The answer is yes. It is contained in the following Theorem. 
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Theorem 1. For any f G L 2 (S 2 ), an arbitrary small e > 0, any fixed a = «o £ S 2 , any 
fixed k = k > 0, and any bounded domain D C M 3 , there exists a (non-unique) potential 
q(x) G L 2 (D), such that (JHJ) holds. 

The relation between the particle distribution density and the potential q is explained 
in Section |2] and Section El In Section 0] we give an algorithm for calculating such a 
potential. Numerical results are presented in Section |S| Our solution of this problem is 
based on our earlier results on wave scattering by small bodies of arbitrary shapes, see 
Ramm (2005b), as well as Ramm (2006a,b). 



2 Scattering by many small particles 

If many small particles D m , 1 < m < M, are embedded in D, u \s m — 0, where S m is the 
boundary of D m , then the scattering problem is: 

M 

[V 2 + k 2 - q (x)]u = in R 3 \ (J D m , (9) 

m=l 

M | 5m =0, m = l,...,M, (10) 

e ikr ( 1 \ , , x 



r 









r = \x\ 







u(x) = u Q (x) + A(j3, a) ho-, r = \x\ — > oo, j3 = — , (11 



and the solution u(x) is called the scattering solution. Here no(x) is a given refraction 
coefficient, n (x) > in M 3 , u is the acoustic pressure, 

qo(x) := k 2 [l — no(x)] = in D' . 

Under the assumptions of Section d we can have d A. We also assume that the 
quantity a/d 3 has a finite non-zero limit as M — > oo and a/d — > 0. More precisely, if 
C m is the electrical capacitance of the conductor with the shape D m , then we assume 
the existence of a limiting density C(x) of the capacitance per unit volume around every 
point x G D: 

C(x)dx, (12) 



lim V" C m = / * 



M 

D m cD 



where D C D is an arbitrary subdomain of D. Note that the density of the volume of 
the small particles per unit volume is 



a 3 \ a 



One can prove (see Ramm (2005b), p. 103) that in the limit M — ► oo the function u solves 
the equation 

[V 2 + k 2 — q(x)]u = in R 3 , q(x) = q (x) + C(x), (13) 
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where C(x) is defined in f)12|). and A(ft,a) in (|TT|) corresponds to the potential q(x) (see 
also Marchenko and Kruslov (1974), where similar homogenization-type problems are 
discussed) . 

If all the small particles are identical, Co is the capacitance of a conductor in the 
shape of a particle, and N(x) is the number of small particles per unit volume around 
point x, then, up to the quantity of higher order of smallness as | — - > 0, we have: 

C(x) = N(x)C . 

Therefore, 

q(x) = q (x) + N(x)C , 

and 

N(x) = 9(x) - go(x) . (14) 
Co 

Thus, one has an explicit one-to-one correspondence between q(x) and the density N(x) 
of the embedded particles per unit volume. 

Remark 1. If the boundary condition on S m is of impedance type: 

u N = (u on S m , 

where N is the exterior unit normal to the boundary S m , and ( is a complex constant, 
the impedance, then the capacitance Cq in formula (9) should be replaced by 

1 + C\s\ 

where \S\ is the surface area of S, and the corresponding potential q(x) will be complex- 
valued, see Ramm (2005b), p. 97. 



3 Scattering solutions 

To establish Theorem ^ recall that for a fixed k > the scattering problem (JT}-© is 
equivalent to the Schrodinger scattering problem for the potential q(x): 

f £ ik\x— y\ 

u q = u - g(x,y)q{y)u q {y)dy, g(x,y) := —. (15) 

Jd 4?r F - y\ 

for which the scattering solution u = u q is the unique solution. The corresponding 

scattering amplitude is 

A(a',a) = — — / e~ ika ' x q(x)u q (x,a)dx, (16) 
4?r Jd 

where the dependence on k is dropped since k > is fixed. 
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If q is known, then A := A q is known. Let q G L 2 (D) be a potential and A q (a', a) be 
the corresponding scattering amplitude. Fix a G S 2 and denote 

A(J3) := A q (a',a), a' = f3. (17) 

Then 

A(3) = -— [ e- ik(i - x h(x)dx, h(x) := q(x)u„(x,a). (18) 
4tt J d 

Our goal in Theorem [T]is to find a potential q for which © is satisfied. First, we find 
an h(x) that satisfies 

\\f((3) + ^J D e-^ x h(x)dx\\ L 2 (s , ) <e. (19) 

The existence of such an h follows from the following Theorem. 
Theorem 2. Let /(/?) G L 2 (S 2 ) be arbitrary. Then 



inf 

h&L 2 {D) 



f(P) + — f e' ik/3 - x h(x)dx 
4vr J D 



= 0. (20) 

L 2 (S 2 ) 



Proof of Theorem^ If ()2L)jl fails, then there is a function f{(3) G L 2 (S 2 ), f ^ 0, such 
that 

/ / e" i/c/3 ^(x)rfx = VheL 2 (D). (21) 

is 2 in 



This implies 



<^>(x) := / dpf(P)e~ ikl3 - x = VxGD. (22) 
is 12 

The function <£>(x) is an entire function of x. Therefore (J22j) implies 

<p(x) = ViG M 3 . (23) 

This and the injectivity of the Fourier transform imply f(0) = 0. Note that (p(x) is the 
Fourier transform of the distribution f(/3)5(k — X)X~ 2 , where 5(k — A) is the delta-function 
and \f3 is the Fourier transform variable. The injectivity of the Fourier transform implies 
f({3)\- 2 5(k - A) = 0, so f(J3) = 0. Theorem |2 is proved. □ 

To find an h that satisfies (fEJJ) one can proceed as follows. Let {Y^(/9)}^ , Yi = Yg >m , 
— I < m < £, be the orthonormal in L 2 (S 2 ) spherical harmonics, 



Y e>m (-(3) = (-l)% m ((3), Y t>m (P) = (-iy +m Y itm (P), (24) 



h(r) := [^) 1/2 J i+h (r), (25) 
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where J? are the Bessel functions and the overbar stands for the complex conjugate. It 
is known that 

e -ik0. x= J- An{-iYh{kr)Y^)Y^ m {(5), x° := ^. (26) 
e=o,-e<m<e ' ' 

Let us expand / into the Fourier series with respect to spherical harmonics: 

f(P)= E ft,mY e ,m(P)- (27) 

l=Q-l<m<t 

Choose L = L(e) such that 

ElA-l 2 ^ 2 - ( 28 ) 

e>L 

With so fixed L, take /?^ m (r), < £ < L, — £ < m < £, such that 

i ji 

km = -{-if{^) J r 3 / 2 J eH (kr)h itm (r)dr, (29) 

where b > 0, the origin O is inside D, the ball centered at the origin and of radius b 
belongs to D, and he tTn (r) = for r > b. There are many choices of /i£ >m (r) which satisfy 
(J2SJ). If (|25jl and (f29~J) hold, then the norm on the left-hand side of (j2Ti|) is smaller than e. 
A possible explicit choice of /i£, m (r) is 

f -(-iY 7 =M - - £ < L, 
he, m =\ V^W*> " (30) 

[0 , £ > L 

where g^ u {k) := J x^ + ^ J v {kx)dx. This integral can be calculated analytically, see 
Bateman and Erdelyi (1954), formula 8.5.8. We have assumed that h(x) = for \x\ > 1, 
and b = 1 in (J2HJ) . Finally, let 

/i(x) = ^J/i| im (r)y^ )m (a'). (31) 

1=0 

This function satisfies inequality (19) by the construction. 

4 Reconstruction of the potential 

In the previous section we have shown how to find a function h e L 2 (D) that satisfies 
(JTHJl . In this section a potential g satisfying the conditions of Theorem ^ is constructed 
from such an h. The possibility of such a reconstruction follows from the following result. 
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Theorem 3. Let h G L 2 (D) be arbitrary. Then 

inf \\h — qu q (x, a)\\ = 0. (32) 

q&L 2 (D) 

Here a G S 2 and k > are arbitrary, fixed. Moreover, if \ \h\\ L 2^ is sufficiently small, 
then there exists a potential q such that 

h(x) — q(x)u q (x,a). (33) 

This Theorem follows from Lemma^and Lemma Instated and proved below. For con- 
venience let us summarize the method for finding a potential q satisfying the conditions 
of Theorem ^ 

Method for Potential Reconstruction. 

Let e > 0. 

Step 1. Given an arbitrary function /(/?) G L 2 (S 2 ) find h G L 2 (D) such that (fiUj) holds. 
This can be done using (J3DJ). Let 

L 

h(x) = Y J Km{r)Y^ m {a') 1 (34) 

1=0 

where L = L(e), see (30). 

Step 2. Use h, obtained in Step 1, to find a potential q G L 2 (D) satisfying 

\\h — qu q (x, a) || < £. 

For / with a sufficiently small norm ||/(^3)||i^(s2) such a potential q can be found using 
formula (|3*7)l. see below. Formula (|3*7|) can be used for any / for which condition (|3T)|) 
holds. 

Step 3. This potential g generates the scattering amplitude A{0) at fixed a and k, such 
that 

||/(^-^)|U 2(s2) <C£ 

holds for some constant C, independent of e. 
Indeed, let || • || := || • ||x^(S 2 )- Then 

H/C9) - A q {(3)\\ = 11/09) + ^ / e- ifc ^g(x) M (x)^|| 



D 

e\D\ 



47T 

This concludes the proof of Theorem ^ 



<\\f(P) + ^ / e^/MfaH + 

e|D| . . 
< e H — - — , LD =meas-D. 



47T 



(35) 
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Lemma 1. Assume that sup x6D 



J D ghdy 



< 1, or, more generally, that 



sup 

rrG-D 



uo(x) - / g(x,y)h(y)dy 



D 



> 0, g = g(x,y):-- 



Then equation (}3*3*j) has a unique solution: 

h(x) 



q(x) 



u (x) = e 



ika-x 



e ik\x-y\ 

4:7i\x — y\ 



q E L\D) 



(36) 



(37) 



Remark 2. It follows from Theorem® and the discussion afterward that f and h are 
proportional, so that if \\f\\ L 2^ S 2^ is sufficiently small, then \\h\\tfm\ is small, and then 
condition (}36|) is satisfied. 

Proof of Lemma QJ The scattering solution corresponding to a potential q solves the 
equation 

u = u - g(x,y)q(y)u(y)dy, u := e ika ' x . (38) 



D 



If h(x) = q(x)u q (x, a) holds, i.e., if h corresponds to a q e L 2 (D), then u = u — f D ghdy. 
Multiply this equation by q and get 

q(x)u(x) = q(x)u (x) - q(x) / g(x,y)h(y)dy. 

Jd 

Using (|3~3~)) and solving for q, one gets (pTTf). provided that (J3H)) holds. Condition (|3lfl) 
holds if ||/i|| L 2( D ) is sufficiently small. One has 



g(x,y)hdy 



D 



and 



dy 



4tt V^d la - y\ 
where a = 0.5diamD. If, for example, 



\x - y\ 



< 



\h\ 



L 2 (D) 



L 2 (D), 



(39) 



47T 



\Hl 2 (D) < i, 



then condition holds, and formula (|37|) yields the corresponding potential. This 
explains the role of the "smallness" assumption. □ 

Remark 3. If ()36j) fail, then formula (J37|) may yzeZd cg^ L 2 (D). As long as formula 
(|37|) yields a potential q G L P (D), p > 1, otzr arguments essentially remain valid. In 
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our presentation we have used p = 2 because the numerical minimization in L 2 -norm is 
simpler. 

The difficulty arises when formula (|37jl yields a potential which is not locally inte- 
grate. Numerical experiments showed that this case did not occur in practice in several 
test examples in which the "smallness" condition was not satisfied. 

We prove that a suitable small perturbation h$ of h in L 2 (D)-norm yields by formula 
(J37|) a bounded potential q$. This means that the "smallness" restriction on the norm of 
f is not essential. 

Lemma 2. Assume that h is analytic in D and bounded in the closure of D. There exists 
a small perturbation hg of h, \\h — hsW^fD) < such that the function 



Outline of proof. Suppose that for a given h G L 2 (D) condition is not satisfied. 
Let us approximate h by an analytic function hi in D, for example, by a polynomial, so 
that 



Denoting hi by h again, we may assume that h is analytic in D and in a domain which 
contains D. We prove that it is possible to perturb h slightly so that for the perturbed 
h, denoted kg, condition is satisfied, and formula (J37)l yields a potential q$ G L 2 (D), 
for which inequality (jHJ) holds, see Ramm (2006c) for details. 

Finally we make some remarks about ill-posedness of our algorithm for finding q 
given /. This problem is ill-posed because an arbitrary / G L 2 (S 2 ) cannot be the 
scattering amplitude A q (/3) corresponding to a compactly supported potential q. Indeed, 
it is proved in Ramm (1992), Ramm (2002), that A(f3) is infinitely differentiable on S 2 
and is a restriction to S 2 of a function analytic on the algebraic variety in C 3 , defined by 
the equation (3 ■ (3 = k 2 . Finding h satisfying (fT§j) is an ill-posed problem if e is small. It 
is similar to solving the first-kind Fredholm integral equation 



whose kernel is infinitely smooth. Our solution (}30|) shows the ill-posedness of the problem 
because the denominator in (JHOJ) tends to zero as I grows. Methods for stable solutions 
of ill-posed problems (see Ramm (2005a)) should be applied to finding h. If h is found, 
then q is found by formula (pT7)l. provided that (J5fij) holds. If does not hold, one 
perturbs slightly h according to Lemma 121 and get a potential q$ by formula (|37j) with 
h$ in place of h. 



h s (x) 



qs ■= 



M x ) ~ J D d(x,y)h 5 (y)dy 



is bounded. 
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5 Numerical results 



In this section we present results of numerical experiments for a design of the material 
capable of focusing the incoming plane wave into a desired solid angle. First, let us note 
that a direct implementation of the algorithm presented in the previous sections produces 
potentials q with large magnitudes (of the order of 10 5 ) in our examples. This happens 
because of the ill-posedness of the inverse scattering problem. To remedy this situation 
we have introduced an additional step in the potential reconstruction algorithm stated 
in Section HJ 

Step lb. Let h^ m be the coefficients of h obtained according to (|30jl . Bound the magni- 
tudes of the coefficients by a predetermined constant T > 0, that is 

hi )T n if \hi >m \ < T 

if \h hm \>T. 

Let 

L 

h(x) = Y,K,m(r)Y e , m (a'). (40) 

1=0 

The bound on the function h has the effect of bounding the potential q. This procedure 
regularizes the ill-posedness of the reconstruction process as discussed at the end of 
Section |3J The ill-posedness manifests itself in the divergence of the series (|34|) with 
L = oo, when the regularization we have used by introducing the bounding constant, is 
not applied. However, from the numerical observations, the series (I40|) practically did 
not change with the increase in L for L > 5. As expected, an increase in the value of T 
improves the precision of the approximation of the desired scattering amplitude /, but 
it also increases the magnitude of the potential q. A reduction in the value of T leads to 
a deteriorating approximation. 

In all the experiments the incident direction a = (0,0, 1), k = 1.0 and L = 6. The 
domain D is the ball of radius 1 centered in the origin. In our first numerical experiment 
the goal was to focus the incoming plane wave into the solid angle < 6 < tt/4, where 9 
is the polar angle measured from the incident direction a = (0, 0, 1). Figure 1 shows the 
cross-section through the incident direction of the desired (dotted line) and the attained 
absolute value (solid line) of the scattering amplitude. 

Figure 2 shows the contour plot of the absolute value of the recovered potential q in 
a cross-section through the z-axis. The darker colors correspond to the larger values of 
\q\. In this experiment the maximum of the absolute value of the potential q was about 
1290 corresponding to the bounding constant T = 100. This value of T was found by 
examining numerical results with larger and smaller values for this constant. For smaller 
T the resulting radiation pattern has smaller magnitudes, i.e. the plot of its absolute 
value is located closer to the origin. For larger values of T the maximal value \q\ of the 
potential approaches the order of 10 5 . 



10 



Figure 1: Attained (solid line), and targeted (dotted line) scattering amplitude /(/3) in 
experiment 1. 

Similarly, Figures 3 and 4 show the results of the numerical experiment aimed at 
focusing the same incident plane wave into the solid angle 0.2ti < 9 < 0.57T. The 
maximum of |g| was about 1840 in this case, corresponding to the bounding constant 
T = 800. This value of T was found experimentally as above. For smaller values of T 
the resulting radiation pattern has a significant component in the region \9\ < 0.2n, i.e. 
it produces a poor approximation for the desired scattering amplitude. 

6 Conclusions 

A method is developed for finding the number N(x) of small acoustically soft particles 
to be embedded per unit volume around every point x in a bounded domain, filled with 
a known material, in order that the resulting new material has the desired radiation 
pattern. Any wave field, not necessarily acoustic wave field, which satifies equations 
(9)-(ll) is covered by our theory. 

On the boundary of each acoustically soft particle the Dirichlet condition holds. 

The method is justified theoretically. Numerical examples of its application are pre- 
sented. The ill-posedness of our problem is discussed and a regularization method for its 
stable solution is proposed and successfully tested numerically. 

The direct application of the derived formula (J30|) may lead to large values of q. To 
remedy this situation the coefficients h^ m are bounded. This is a way to handle the 
ill-posedness of the inverse problem. The resulting algorithm exhibits a stable behavior. 
It serves as a regularizing algorithm for solving the original ill-posed problem. Numerical 
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Figure 2: Contour plot of the potential q in experiment 1. 
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Figure 3: Attained (solid line), and targeted (dotted line) scattering amplitude f(/3) in 
experiment 2. 
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Figure 4: Contour plot of the potential q in experiment 2. 
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results show that the method can produce materials with the desired focusing properties 
under the limitation that the desired radiation pattern /(/?) is well approximated by a 
short series of spherical harmonics. 
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